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Abstract 

CL| We consider some aspects of nonlocal modified gravity, where nonlocality is of the 



type RT{\3)R. In particular, using ansatz of the form DR = cR^ , we find a few R{t) 
solutions for the spatially flat FLRW metric. There are singular and nonsingular bounce 
solutions. For late cosmic time, scalar curvature R{t) is in low regime and scale factor 
a{t) is decelerated. R{t) — satisfies all equations when fc = — 1. 

1 Introduction 

General theory of relativity was founded by Einstein at the end of 1915 and has been 
successfully verified as modern theory of gravity for the Solar System. It is done by the 
Einstein equations of motion for gravitational field: i?^^ — ■^Rgf^i, — KT^,y, which can be 
derived from the Einstein-Hilbert action S = J y/^Rd^x + / 

Attempts to modify general relativity started already at its early times and it was 
mainly motivated by research of possible mathematical generalizations. Recently there 
has been an intensive activity in gravity modification, motivated by discovery of ac- 
celerating expansion of the Universe, which has not yet generally accepted theoretical 
explanation. If general relativity is theory of gravity for the Universe as a whole then it 
has to be some new kind of matter with negative pressure, dubbed dark energy^ which 
is responsible for acceleration. However, general relativity has not been verified at the 
cosmic scale (low curvature regime) and dark energy has not been directly detected. This 
situation has motivated a new interest in modification of general relativity, which should 
be some kind of its generalization (for a recent review of various approaches, see [T], 
and for renormalizability [2]). However there is not a unique way how to modify general 
relativity. Among many approaches there are two of them, which have been much inves- 
tigated: 1) f{R) theories of gravity (for a review, see [3]) and 2) nonlocal gravities (see, 
e.g. [Ills] and references therein). 

In the case of f{R) gravity, the Ricci scalar R in the action is replaced by a function 
f{R). This is extensively investigated for the various forms of function f{R). We have 
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had some investigation when f{R) — i? cosh "^^^^ and, after completion of research, the 
results will be presented elsewhere. 

In the sequel we shall consider some aspects of nonlocal gravity. Nonlocality means 
that Lagrangian contains an infinite number of space-time derivatives, i.e. derivatives up 
to an infinitive order in the form of d'Alembert operator □. In string theory nonlocality 
emerges as a consequence of extendedness of strings. Since string theory contains gravity 
as well as other kinds of interaction and matter, it is natural to expect nonlocality not only 
in the matter sector but also in geometrical sector of gravity. On some developments in 
cosmology with nonlocality in the matter sector one can see, e.g., [6l[7l[8j[9] and references 
therein. In the next section we shall discuss a nonlocal modification of only geometry 
sector of gravity and its corresponding cosmological solutions (on nonlocality in both 
sectors, see [lO]). 

2 On a Nonlocal Modification of Gravity 

Under nonlocal modification of gravity we understand replacement of the Ricci curvature 
R in the action by a suitable function F{R, □), where □ = -^^^dfj,^J—gg^'^dy. 

Inspired by [5 (for recent developments, see [HI [12]), we consider nonlocal Lagrangian 
without matter in the form 

oo 

which was proposed in [1^, where J^(n) = /nQ" and c is a constant. By variation of 

n=0 

the Lagrangian ([!]) with respect to metric g^^ one obtains the equation of motion for <7^y 

+00 1 
n=l 1=0 

+ d^U^Rd^.W-^-^R - g^^igP^dpU^RdaW'-^-^R + U^RH'-^R)) 

- A-KGg^^,cRF{U)R + IQ-KGciD^d,, - g^^n)F{U)R. (2) 

The trace of ^ is also a useful formula and it is 

-t-C30 n— 1 

y /„ V(aun'i?5^n""^"'i? + 2n'i?n"-'i?) + &ut{u)r = —-. (3) 

^-^ ^-^ SvrGc 

n=l 1=0 

We mainly use the spatially flat (homogeneous and isotropic) Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric ds^ — —dt^ + a^{t){dx'^ + dy"^ + dz^). Investigation 
of ^ and finding its general solution is a very difficult task. Hence it is important to 
find some special solutions. To this end some ansatze of the form OR = cK^ seem to be 
useful. In the sequel we construct a few such ansatze. 
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2.1 Case OR = rR 

At the beginning, to illustrate method, we investigate ansatz of the simplest form: OR = 
rR. For this ansatz, where r is a constant, we have 

□"i? = r"R, T{a)R = T{r)R. (4) 
In the FLRW metric □ = —d^ — 3Hdt and ansatz OR — rR becomes 

R + 3HR + rR^O, (5) 
where = ^ is the Hubble parameter. Replacing 

R^6{H + 2H^) (6) 

in Eq. ([5]) we get 

H + 4:H^ + 7HH + 12H'^H + r{H + 2H'^) = 0. (7) 



A solution of this equation is 



This implies scale factor a{t) — C2\/\'2t + Ci\ and acceleration a = ^, , 

\2t+Ci\y/ \ 2t+Ci\ 

where C2 > 0, Ci € R. Calculation of R by expression (jS]) gives R — and it is consistent 
with other formula containing R, including 

It is natural to take constant Ci = 0, because it yields symmetrical solutions with 
respect to t = 0. Result a{t) = C2\/\2t\ is an example of the symmetric singular bounce 
solution. 

Note that the above solutions hold also when r = in the ansatz, i.e. OR ~ 0. More 
general ansatz OR = rR + s was considered in [S] . 

2.2 Case DR = qR^ 

The corresponding differential equation for the Hubble parameter is 

H + 4ij2 + 7HH + 12H^H + 6g(ij^ + m'^H + AH^) = (9) 

with solution 

' 3 t + Ci (27;+l)(4r;-l)' ' ^ ' 

Another solution is H = ^ t^Ci ^i^h arbitrary coefficient q, what is equivalent to the 
ansatz OR = rR with R — 0. 

The corresponding scalar curvature is given by 

2 (27?+l)(4ry-l) 
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It is interesting that □"i?,! = when n e N. This can be shown by mathematical 
induction by the following way. It is evident that DRi = 0. Suppose that — 0, 

then □"+ii?„+i = □□"i?„ + i6"+i0 n»Qj^^ ^ 0. 

This □"i?„ = property simplifies the equations considerably. For this special case 
of solutions trace equation ([3]) effectively becomes 

V fk V(5.n'i?a^n'^-i-'i? + 2n'i?n'^-'i?) + 6nj"(n)i? — — , (12) 

^-^ ^-^ SttGc 

k=l 1=0 

where 

n-l 

J^{D)R = ^ hW'R. (13) 

fc=0 

In particular case n — 2 the trace formula becomes 

- foR^ + hi-R^ + —R^) + f2{-—RR^ + -^R^) + h{-^R^R^) 
35'' ■' ^ 35 '^-'^"^ 35 ^ 1225 ^^•''^^ 1225 ' 



R 



(14) 



SvrGc' 

2.3 Case = c„i?""+^ 

We considei0 another ansatz of the form □"i? = c„i?""''"''; where a and (3 are constants, 
and n € N. From the equalities 

= c„((an + /3)i?""+'^-ini? - (an + /3)(an + /3 - 

= c„(an + /3)(cii?""+"+2^-i - (an + /3 - ^ c„+ii?""+"+'^ (15) 

we get the following conditions: 

an + a + 2/3 - 1 = an + a + /3, (16) 
7?2 (17) 
c„+i = c„(an + /?)(ci — an — /3 + 1). (18) 



Equation (fT6|) implies that /3 is equal to 1. Sequence c„ is defined by (fT8|) and can be 
explicitly written (/3 = 1) as 



ri-l 



= ci (a/s + l)(ci — a/c), (19) 



where ci is a constant. General solution of equation ([TT]) is of the form 

i?(<) = 22/"(a(±<-di))"'/", rfiGM (20) 



^B. D. thanks A. S. Koslielev for suggestion of ansatz ~ _R"+^. 
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with arbitrary constant di . 

Case a = 1. In the case a = 1 the coefficients c„ are given by c„ = where ci 

is the first element. Putting a = 1 into equation ((20|) one obtains 



{t-di) 

The corresponding expressions for H{t) and a{t) are: 



= (21) 



Hit) ^ i^ + Vm<^At-d.)^-~V^ + ^^ (22) 
i2{t-di) (d2 (^-c^l)^/^ + l 



3-n/57 

ait) = ds {t ^ di) 1^ Vc^2(i-di)^ ^ +1, (23) 

where di, d2, rfa are arbitrary real constants. 

The function ait) has a vertical asymptote at the point t = di. If di < then a(t) > 
for all t > 0. For large values oft, a(t) is asymptotically equivalent to t 2\/^+T2 (3-^/57) 

^0.879 

d3rA(-21- V57) _ 5-) ^2y2yf _ 48rf2rV¥ - V57 - 5 

24 (d2^^/T + 1) 



°(^) = ^ : 7172 (24) 



where T = t — di. The expansion is accelerated for d2it — (ii)V^ < ^^21 \- and 

v57— 5 v57— 5 

it is decelerated otherwise. 

Note that this ansatz □"i? = CnR"^^ for n = 1 coincides with ansatz DR = q^i?^, 

when 77 = — because then one can take ci = q,, = In this particular 

case they have the same scalar curvature R and the same Hubble parameter for ^2 = 0. 
However, apart from this special case 77 = constant q,, is different of ci = -9±%/57 _ 
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Figure 2: Scale factor a{t) given by (p7|) ioi di = d2 = 2 and 



da - jQ. 



Case a 



From 



Putting a — ^ into equation (|20p we obtain 



256 



(t + di)' 



we obtain 



32 / 32 ~ 

-512\/3d2e-^<''i+') + 3(i + rfi) f 32^26 ^(<^i+*) + %/3 



and then 



a{t) = dae v^<'i+') y/di + t\J 32^26 + Vs, 



where di, (i2, c?3 are some real constants. 
The corresponding acceleration is 



a{t) 



12 (di + 7/2 (32^26 -^(''1+*) + \/3) 3/2 
(^-6dit - 3dl + 32V3di - 3i^ + 32^3^ + 256 
3 (edit + 3dl + 32V3di + 3<2 + 32V3t - 256 
192\/3d2e-^<'i+*' (di + i - 16) {di + t + 16) ) 



1024d^ ev^<<=i+') 



(25) 



(26) 



(27) 



(28) 



The acceleration is positive for t < to and negative for t > to-, where to is the zero of 
a{t). For large values of t, d{t) converges to 0. 
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Case a — 2. In the case a = 2 only one integration can be performed and it gives 



Hit) 



V3 ; "'^^ V V3 



^ -y/F^i^o (^) + ^/3i.: (^) - yi^i^. (^) ^^^^ 

Ji and are modified Bessel functions of the first and the second kind, respectively, and 
di, d2 are real constants. 

Case a = —2. For a — —2 we obtain expression for Hit) involving Airy functions 

{Pi (daAi' ( {r~\ (t -d,))+ Bi' ( (t -di))) 
H{t) = -^-4^ . . (30) 



(d2Ai(^(t-rfi)) +Bi(^(i-di))) 



3 Concluding Remarks 

In this article we presented three ansatze, two of them are quite new and one can be 
adjusted so that □"i? = 0. These two ansatze have solutions for scalar curvature of the 
form R = (^f^^_^j2 1 which satisfy all but extended Einstein equations ^ and related trace 
formula It is a consequence of the quadratic form in R of the Lagrangian ([T]). However 
these ansatze are promising for some new nonlocal Lagrangians, which investigation is in 
progress. 

It is worth mentioning that all the above ansatze contain solution i? = 0, which 
satisfies all (including ([2]) and ([3])) equations with curvature constant fc = — 1. Namely, 
for R= 0, Eq. ([2]) reduces to G^i/ = and it gives 

- = 0, (- +- = 0. 31 

If A: = one has only static solution a — constant. However, when k ~ —1 then a{t) = \t\ 
and it contains a crunch preceding to a big bang. 

Above considered ansatze may be also useful in analysis of some other nonlocal gravity 
and cosmology models. Further investigation of nonlocality governed by the Riemann 
zeta function in p-adic strings dynamics |14| extends interesting cases and can give new 
insights. 
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